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Abstract 

We study the stationary solutions of the Gross-Pitaevskii equation 
that reduce, in the limit of vanishing non-linearity, to the eigenfunc- 
tions of the associated Schrodinger equation. By providing analytical 
and numerical support, we conjecture an existence condition for these 
solutions in terms of the ratio between their proper frequency (chem- 
ical potential) and the corresponding linear eigenvalue. We also give 
approximate expressions for the stationary solutions which become ex- 
act in the opposite limit of strong non-linearity. For one-dimensional 
systems these solutions have the form of a chain of dark or bright 
solitons depending on the sign of the non-linearity. We demonstrate 
that in the case of negative non-linearity (attractive interaction) the 
norm of the solutions is always bounded for dimensions greater than 
one. 

PACS: 03.65.Ge, 03.75.Fi, 47.20.Ky 



1 Introduction 

Recent achievement of Bose-Einstein condensation (BEC) in gases of alkali 
atoms has generated an impressive amount of experimental and theoretical 
works I], . In these systems the condensate is usually described by the so- 
called Gross-Pitaevskii equation (GPE), a Schrodinger equation with a local 
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cubic non-linear term which represents the interaction among the bosons in 
a mean field approximation. GPE effectively reproduces the ground state 
properties of a condensed boson gas confined by an external potential at 
zero temperature ||. In the framework of linear response theory, the mean 
field approximation also allows to evaluate the spectrum of the excitations 
in presence of an external time-dependent perturbation p, |J. On the other 
hand, GPE also appears in the description of other physical systems, like 
nonlinear optics |5[], molecular physics ||, etc.. 

In this paper we study general properties of the stationary solutions of 
GPE both in the case of repulsive and attractive interaction. Besides its 
mathematical interest, this study is relevant in the search of the so-called 
vortex states and, in general, in understanding the dynamical properties of 
condensates. We choose to work in the grand-canonical ensemble, that is 
we fix the chemical potential fi of the system, i.e. the proper frequency 
for the time evolution, and derive the number of particles corresponding to 
each stationary solution. In particular, we study the stationary solutions 
of GPE that have a linear counterpart, in the sense that they reduce to 
the eigenstates of the linear Schrodinger equation which is the limit of GPE 
for vanishing interaction. For these states we conjecture an existence con- 
dition which depend on the ratio between the chemical potential and the 
corresponding eigenvalue of the associated Schrodinger equation. We give 
a proof of this conjecture for the node-less state of a system in presence of 
a general external potential and verify it for the exactly solvable case of a 
one-dimensional, infinitely deep, square well. We also provide numerical ev- 
idence of the validity of the conjecture by studying systems with harmonic 
potentials in different dimensions. 

As a consequence of the above conjecture, we find that in the case of 
attractive condensates there exists a range of the chemical potential \x in 
which the node-less stationary solution does not exist and the lowest-energy 
state has one or more nodes. This may be relevant for the observation of 
stable vortex states. 

We also study the limit of strong non-linearity of GPE obtained for large 
values of the modulus of the chemical potential \x. In this limit a Thomas- 
Fermi approximation holds for repulsive systems, while for attractive systems 
the solutions become independent of the external potential. In the one- 
dimensional case the corresponding approximate solutions have the form of 
a chain of dark or bright solitons depending on the sign of the non-linearity. 
We use these asymptotically exact expressions to establish that the number 
of particles in the ground state of an attractive condensate is always bounded 
for dimensions greater than one in agreement with previous numerical results 



2 The linear limit 

We consider the Gross-Pitaevskii equation || [| describing, in the mean 
field approximation, a system of interacting particles confined by an external 
potential V(x) 

with x G M d . The constant Uq is positive (negative) in the case of repulsive 
(attractive) interaction. Equation (|l|) has two conserved quantities, namely 
the number of particles (squared norm) 



\^(x,t)\ 2 dx 



and the energy 



r ^2 



^-\VV(x,t)\ 2 + ^|*0M)| 4 + V(x)\*(x,t)\ 2 



(2) 



dx. (3) 



The stationary states of Eq. (H), ^(x, t) = exp (— Jr//i) ip(x), where \i is 
the chemical potential are determined by the equation 



2m 



V 2 ij(x) + U \ij(x)\ 2 ij(x) + V(x)i){x) - nif>(x) = 0, 



(4) 



i.e. as critical points of the grand-potential functional 

m = If^iv^p + ^i^r+^^-zi)^)! 2 

= E[il>]-iiN[il>]. 
It is simple to show |l(J that if ip is a solution of (§) then 



da; 



~ I Mx)\ 4 dx. 



(5) 



(6) 



We will look for the solutions of @ corresponding to a given chemical 
potential \i. In this paper, we concentrate on solutions which admit a linear 
counterpart in the sense that they reduce, in a proper limit, to the eigen- 
functions of the associated linear problem 



h 2 

— V 2 (f) n (x) + V(x)(f) n (x) - £ n (pn(x) = 0. 



(7) 
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Here we suppose that So < £\ < • . • < £ n and {<p n (x)} is a hortonormal base 
with 0o (x) positive and bounded. Solutions without linear counterpart will 
be discussed in another paper JiT]. 

By substituting ip(x) = ^ N(fi)x(x) in (|J) with \\x\\ = 1, we have 

-^V 2 xN + U N(ji)\ X (x)\ 2 x(x) + (V(x) - y) X (x) = 0. (8) 

If the number of particles is sufficiently small, the nonlinear term in (§) can 
be neglected and x approximated by n . By substituting x with (fi n , Eq. 
(P) provides the following relation between the chemical potential /i and the 
corresponding norm N(n) 



n~e n + u N(n)\\(i>. 



2 1 1 2 



(9) 



Equation (Q) suggests the following conjecture for the existence of solutions 
of (f|) with linear counterpart 

Conjecture. ForU$ > {Uq < 0), solutions with linear limit ~ ^/N(fi)(j) n 
exist only if ji> £ n (/i < £ n ). Moreover N(n) — > for /i — > £ n . 

In Appendix [A| we give a general proof of this conjecture in the case n = 0. 

The conjecture can be verified analytically in the case of a 1-dimensional 
system confined in a box of size L, i.e. with 



r o 


\x\ 


[ oo 


\x\ 



x\ < L/2 



V{x) 



For this problem the solutions of (|]) are known 
are given by the Jacobi elliptic functions 



(10) 

In the case U > they 



where 



ip n (x) = A sn 2(n + l)K(p) 



K(p) 



L + 2 



=d0 



12) 



'o a/1 — p sin 2 # 

is the complete elliptic integral of the first kind with modulus p G [0, 1], and 



n — 0, 1, 2, By substituting (|Tl|) into one finds the conditions 

ft 2 



A 



mU L 2 

h 2 p + 1 
ml? 2 



p(2(n + l)K(p)) 2 , 
(2(n + l)K») 2 . 



(13) 
(14) 



The number of particles and the energy are given by 

^)~(2(n + l)K(p)) 2 (l-||), (15) 



E(fJ.) = N£ 



7T J 



E(p) 

K( P ) 



(16) 



where 



E(p) 



1 — p sin 2 6* ci^ 



(17) 



is the complete elliptic integral of the second kind with p determined in terms 
of \i by Eq. (fU]). Since K(p) increases monotonously from K(0) = tt/2, for 
a given n Eq. (PHI) has solution only if 



[n + l) 2 7r 2 fi 2 
2mL 2 



H>£ n = 



which complies with the conjecture formulated above. The same conclusion 
can also be reached by using the theorems of [|10| . 

In the linear limit \i — > S n , the solutions ( |TTD reduce to the eigenf unctions 
of the associated Schrodinger equation 



1 



fi— >£ n 



Sill 



T + n ) ( n + 1 ) 7r 



(19) 



In the opposite limit of strong nonlinearity, fi ^> S n , we get from ( |TTD the 
dark soliton solutions 



(x - x k ) ) , 



fc=0 

L L 

2 n+1 



k. 



Similar results are obtained in the case Uq < 0. The solutions of 
now given by 



ip(x) = A cn 2(n + l)K(p) 



x 1 



+ K(p) 



P 



(20) 

(21) 
) are 

(22) 
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with the conditions 

A 2 = -^L* p(2{n + 1)K{P))2 (23) 

Number of particles and the energy become 

N ^=-^L i2in+1 ^ 2 { p - 1+ m)' (25) 



7T 



P 1 + K(p) 



(26) 



where p is determined by Eq. Since (1— 2p)K(p) decreases monotonously 

for p G [0, 1], the n-node solution exists only if \i < S n as conjectured above. 

For the solutions (]22j) have the same limit dH). For -// > £ n , 

we get the bright soliton solutions 

m*) ^D-^f^P^-^). ( 27 ) 

x fc = + (& + -). 28 

2 n + 1 V 2/ v; 

In Fig. we show the behaviour of iV(/z) evaluated according to ( |i3| ) and 
(|25|) for the states n = and n = 1. Note that N(p) — *> for fi — > £ n . The 
single-particle energy, E(/j,)/N(/j,), for the same states is shown in Fig. ^| 

We have verified the above conjecture with numerical and analytical 
methods also in the case of a quadratic potential V(x) = ^m^ 2 Yli=i x h 
with d — 1,2, 3. For example, in the case d — 1 consider the following 
Ansatz for the solutions of @ 

ip n (x) = a n exp (-7^p) H n (jf) , ( 29 ) 

where H n (x) is the Hermite polynomial of degree n and a n , b n are real 
constants. Extremization of the functional Q with respect to a n and b n leads 
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Figure 1: Number of particles N as a function of the chemical potential ji 
for the one dimensional square well. The solid and dashed lines are given by 
Eqs. (|15D and (|25f) , respectively. The two curves correspond to the states 
n = and n = 1. 



to 



2 n+1 n! ( 8 - v/4 + 15(2n + 1) V 

«n = £ " 77 ~, (30) 

(To I5g n 



ti 2 2 + ^4+ 15(2n + 1)V 
m/i 5(2n + 1)77 



(31) 



where 77 = h 2 uj 2 / ji 2 and g n = H n (x) 4 dx. If £/o > 0, the condition 



-00 

a 2 > implies // > (n + |) faj. If £/ < 0, the same condition leads to 
H < (n + ^) hw. Note that in the linear limit /i — > (n + ^)fkv, we have 
A" oc a 2 n — > and fe 2 — > h/mu. 

Analogously, in the case d = 2 consider the Ansatz 

*•(*»■*>=■■ e ^-3 F (~ n ' |m|+1 'fe)) (£f e,ra °- (32) 

where r 2 = x\ + x 2 ,, tan# = x%/xi and F(n,m,r) is the confluent hyperge- 
ometric function ]T3|. The condition a 2 > is now equivalent to // > £ niTO 
if [To > 0, and /i < £ n>m if U < 0, where £ n m = (2n + |m| + 1) hu are the 
eigenvalues of the associated Schrodinger equation. 
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Figure 2: Single-particle energy £7/iV as a function of the chemical potential 
H for the same states of Fig. [l]. The solid and dashed lines are given by Eqs. 
fll6|) and (p6[), respectively. 



Equation (01) has also been solved numerically with a standard relaxation 
algorithm JTJJ]. In Fig. [3] we show the number of particles obtained as a 
function of the chemical potential \x for the states (n,m) = (0,0), (0, 1) and 
(1, 0) in the case of a two-dimensional quadratic potential. The single-particle 
energy for the same states is shown in Fig. || Similar results are obtained 
for d = 1 and d = 3. 

Figures allow us to emphasize a possibly important consequence of the 
above conjecture. In the case Uq < 0, the node-less solution exists only for 
H < £q. Therefore, in the range Sq < fi < Si the state with minimal energy 
is This implies that controlling the chemical potential it is possible to 
obtain a condensate with a node or a vortex in the ground state. 



3 The strongly non linear limit 

The conjecture discussed so far concerned the behaviour of the solutions of 
Eq. in the linear limit. An approximate expression of the solutions of (Q) 
is possible also in the opposite limit — > oo. Let us consider first the case 
Uq > 0. The repulsive interaction tends to delocalize the solutions so that 



the Thomas- Fermi approximation holds [15, 16 1. In this case the gradient 
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Figure 3: Number of particles N as a function of the chemical potential /x for 
a two dimensional quadratic potential. Solid and dashed lines are obtained 
by the solving numerically Eq. (|J) for U > and U < 0, respectively. The 
three curves correspond to the states (n, m) = (0,0), (0, 1) and (1,0). 

5 I ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 
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Figure 4: Single-particle energy E/N as a function of the chemical potential 
fx for the same states of Fig. |3|. 
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term in Eq. (|]) can be neglected and if) is determined by 

U \i)(x)\ 2 i)(x) + V(x)il>(x) - ^{x) = 0. (33) 
Therefore the ground state solution can be approximated as 

Mx) = I f^vWM , > vf) (34) 

[0 [A < V(x) 

In the one-dimensional case, n-node solutions may be approximated by a 
chain of dark solitons 

ip n {x) = ip (x) PJtanh ( ^-^ (x - x k ) J (35) 



k=l 



with Xk to be determined, for instance by extremizing the functional Q. 

In the case of a quadratic potential the number of particles and the energy 
for the state ( |5^ ) are 

d+2 d+2 

»M = 5^A M J^, (36) 



£M = J V^l-_ y (37) 

where A(d) is the volume of the unitary <i-dimensional sphere. From Eq. ( |36| ) 
we see that N diverges for fi — > oo. Similar results are obtained for other 
potentials. 

In the attractive case Uq < 0, the solutions of (El) tend to localize and the 



Thomas- Fermi approximation fails JT5J. In this case, however, for fi — > — oo 



the potential term VV> becomes negligible and Eq. (||) can be approximated 

as 

_|lvV(^) + U \iP(x)\ 2 il)(x) - ^{x) = 0. (38) 

Recently numerical evidence has been provided that the number of par- 
ticles confined in a two dimensional harmonic potential is limited in the case 
of attractive interaction (7J. This fact can be analytically understood from 
(H). With the change 

x = -j^t (39) 

1>(*) = M M)> ( 4 °) 
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Eq. ([38]) can be rewritten in the adimensional form 



l^im-\m\ 2 m+m = o. m 

Note that in the one-dimensional case, n-node solutions may be approximated 
by a chain of bright solitons 




fc=0 



l) fc sech(^^(*-* fc )) 



Mx) = \ TT > f-l) fe sech v - ^ {x - x k ) (42) 



with Xk to be determined, for instance by extremizing the functional Q. The 
number of particles corresponding to a solution of fl4"I|) is given by 

n^) = [ mx)fd X = ^— * [ \m\ 2 dt 

J U (_ m/i )2 J 

h " >|^r 2 (d), (43) 



m*\Ut 



o 



where Tk(d) = J \(fi($,)\ k d$, is a numerical constant. Therefore we have 

oo d = 1 

lim Nfa)=< ^ r ^(2) d = 2 . (44) 



m|!7o| 

d > 3 



In Fig. H] we show the behaviour of N(fi) in the ground state obtained 
by solving numerically Eq. (^) with a harmonic potential in the cases d = 
1,2,3. We have chosen the following realistic values for the parameters: 
m = 3.818 x 10 26 Kg, lu = 10.0 Hz and, for d = 3, U = Aith 2 a s /m with 
a s = 2.75 x 10~ 9 m fllTf. For d = 2 and d = 1 we set C/ = Anh 2 a s /mL 
and Z7 Q = A-ah 2 a s /mL 2 with L = 10" 5 m and L 2 = 9 x 10" 10 m 2 . The 
numerical results compare very well with the analytical approximations ([36f) 
for f/o > and (|43f) for f/ < 0, respectively. In the case of Eq. (|43|), T 2 (d) 
has been evaluated numerically. We have T 2 (l) = 2.82842, T 2 (2) = 5.85044 
and T 2 (3) = 6.68118. Note that for d = 3, N(fi) has a maximum and vanishes 
for both fi — > — oo and /i — > 3/2hu. This implies that the function fi(N) is 
not single- valued but has two branches in agreement with []18| . 

In Fig. ^] we show the single-particle energy evaluated numerically in the 
same cases of Fig. |^. For fi — > oo the energy diverges for any value of d 
according to the limiting expression fl3"7p. For /i — > — oo the behaviour of 
E(fi) is well described by 

= w " - IWi) (45) 



11 



10' 



10" 



10= 



a. 



jo 



10 J 



10' 



10 



10" 




-30 -20 -10 10 20 



30 



Figure 5: Number of particles N as a function of the chemical potential \i 
in the ground state of a quadratic potential with d = 1 (dotted line), d = 2 
(dashed line) and d = 3 (solid line). The dots + and x are the analytical 
results (BB) and (f43|), respectively. 



which easily stems from Eq. fl4"T|). For d = 2 we have T4 = 2r2 and hence 
-£?(//) vanishes for \i - 



-oo. 



A Existence of a node-less state 

In the following we suppose that the external potential is bounded from below 
and, for simplicity, we take V(x) > 0. We will prove that, with U > 0, a 
solution of (HI) exists if and only if /j, > Sq. The proof of the necessary 
condition is based on the property (H). Let us define the functional 



Qo[ 



^|V#s)| 2 + (V(x) - AMx) 



dx. 



(46) 



We have n[ip] = Q [ip] + \U Q J \^(x)\ 4 dx. If Q [^\ > 0, then ip(x) cannot 
be a solution of (^). The linear problem Q' [(p n ; x] = k n (j) n (x), where 



Q'o[<t>n,x\ 



5Qo[(pn] 



[X 



A V 2 0n( : r) + (V(x) - fi) <j> n {x), (47) 
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Figure 6: Single-particle energy E/N as a function of fi in the ground state 
of a quadratic potential with d = 1 (dotted line), d = 2 (dashed line) and 
d = 3 (solid line). The dots + and x are the analytical results (|37|) and fl45[), 
respectively. 



has the same eigenfunctions of (0) and the eigenvalues are k n = £ n — /i. By 
decomposing a generic ^>(£c) as ^(aj) = Y^Lo c n<fin(x), we obtain 

Qo[1>] = (Q'aW;xU(x)) 




(48 



n=0 



Therefore, if /x < £o we have fco > and Qo[if)] > 0. 

The sufficient condition can be proved with the help of general theorems 
on elliptic differential equations First we look for upper and lower 



solutions of (f|). An upper solution ip u (x) is defined by 

-tJ-V 2 ^) + U Q \^ u {x)\ 2 ^ u {x) + (V(x) - ^ u {x) > 0. (49) 
2m 

For a lower solution ipi(x) the inequality is reversed. If a couple of ordered 
upper and lower solutions exist, i.e. ij) u > fy, then the existence of, at least, 
one solution ip(x) with < tj) < ij) u is guaranteed |B|. It is simple to check 
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that an upper solution is ip u (x) = yJjJUo- As a lower solution we choose 
ipi(x) = e<j) (x) with 



e < mm 



li - S 



U max x |0o(aO| 2 ' v^maxx \(f>o{x) 



(50) 



which ensures that ipi < ip u . 

In the case U < 0, it is possible to prove that a positive solution of (£|) 
does not exist if /i > 8 . Multiplying (£|) by 4>q(x) and integrating, we have 







Therefore, 



ip{x) 



-^-V 2 + Uo\i,{x)\ 2 + {V{x)-^ 
2m 

-7^ 2 + u mx)\ 2 + (v(x)-fi) 



ip(x) dx 
j)o(x) dx 



^(x) [U \^j{x)\ 2 + {S - //)] <t> Q {x) dx. 



(51) 



Uo I (po(x)\ip(x)\ 2 ip(x)dx = (/i - So) / <po{x)i){x)dx 



(52) 



If tp(x) is a positive function, both integrals in (|52|) are positive and for 
fi > S the equality is impossible. 
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